Strength of the interactions in YBa 2 Cu 3 OQ,7 obtained from the inelastic 

neutron-scattering measurements 



Z. G. Koinov 

Department of Physics and Astronomy, University of Texas at San Antonio, San Antonio, TX 78249, USJ^ 

It is widely accepted that: (i) the angle-resolved photoemission spectroscopy (ARPES) data 
produce evidences for the opening of a d-wave pairing gap in cuprates compounds described at 
low energies and temperatures by a BCS theory, and (ii) the basic pairing mechanism arises from 
the antiferromagnetic exchange correlations, but the charge fluctuations associated with double 
occupancy of a site also play an essential role in doped systems. The simplest model that is consistent 
' with the two statements is the t-U-V-J model. We have shown that the inelastic neutron scattering 

| data on YBa 2 Cu 3 Oe.7 [Phys. Rev. Lett. 83, 608 (1999)] combined with the corresponding ARPES 

data allow us to obtain the strength of the on-site repulsive interaction U, as well as the strengths of 
the spin-independent attractive interaction V and the spin-dependent antiferromagnetic interaction 
J. 
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. PACS numbers: 71.10.Fd, 71.35.-y, 05.30.Fk 
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5^ ■ The magnetic susceptibility probed by neutron scattering in cuprates compounds has provided a variety of exper- 
Oh| imental peaks associated with incommensurate and commensurate structure. This neutron data remains of central 
importance in the field of high-Tc superconductivity because the technique could be used to cover momentum and 
frequency range which is wider than that of any other spectroscopies. The neutron data in cuprates compounds are 
. characterized by a very strong dependence on energy, momentum transfer, temperature, and doping. For example, 
the existence of a commensurate peak at ~ 40 meV and incommensurate peaks at ~ 24 meV and ~ 32 meV has 
I , been reported in undcrdoped YBa2CusOe.T^. The simplest and most transparent hypothesis put forward by many 
authors^ is that the commensurate resonance and the incommensurate peaks in cuprates compounds have a common 
origin. This approach (also known as the Fermi-liquid approach) is based on the assumption that the electrons near 
the Fermi surface are strongly correlated due to the on-site Hubbard and nearest-neighbor antiferromagnetic inter- 
actions. Alternatively, according to the stripes model (see Ref. [H and the references therein) the incommensurate 
peaks are the natural descendants of the stripes, which are complex patterns formed by electrons confined to separate 
linear regions in the crystal. 

In this paper we analyze the position of commensurate and incommensurate peaks assuming the Fermi-liquid-based 
scenario which consists of two steps. The first step is to obtain the tight-binding form of mean-field quasiparticle 
energy and the corresponding chemical potential by matching the shape of the Fermi surface measured by the ARPES. 
y—( • Another input parameter for the d-wave superconductivity is the maximum gap which usually is assumed to be equal 
to the ARPES antinodal gap. During the second step of our approach we use the Bcthe-Salpeter (BS) equations to 
identify effective interaction strengths U, V and J that are consistent with both ARPES measurements and the energy 
and the position of the commensurate and incommensurate resonances observed by neutron scattering experiments. 
7— ( . It is known that the YBaCuO is a two-layer material, but most of the peak structures associated with the neutron 
cross section can be captured by one layer band calculations^. The effects due to the two-layer structure can, in 
principle, be incorporated in our approach, but this will make the corresponding numerical calculations much more 
complicated. In the one-layer approximation, the Hamiltonian of the t — J — U — V model contains terms representing 
the hopping of electrons between sites of the lattice, the on-site repulsive interaction between two electrons with 
opposite spins, the attractive interaction (due to phonons) between electrons on different sites of the lattice and the 
spin-dependent Heisenberg near- neighbor interaction (due to short-range antiferromagnetic order): 

H=- Ei,j> tiM,<ri>j,o - A* E*, CT ™i,<r + U n^n iA 

- V T,<i, j>a a' + J E<i,j> Si.Sj, (1) 

where is the chemical potential. The Fermi operator ipj a (V^o-) creates (destroys) a fermion on the lattice site i 
with spin projection a =t,i along a specified direction, and rii j(T — ip\ is the density operator on site i with 
a position vector r^. The symbol Yl<ij> means sum over nearest-neighbor sites. The first term in |T]) is the usual 
kinetic energy term in a tight-binding approximation, where ttj is the single electron hopping integral. We assume 
V > 0, so the fourth term is expected to stabilize the pairing by bringing in a nearest-neighbor attractive interaction. 
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The last term describes the nearest-neighbor spin interaction. The spin operator is defined by Is i = ipj ^ aa'^>i,a> /2, 
where H = 1 and a is the vector formed by the Pauli spin matrices (<j x , a y ,a z ). The lattice spacing is assumed to be 
a = 1 and the total number of sites is N. 

The magnetic neutron scattering directly measures the imaginary part of the generalized spin susceptibility 5 
Xr/(Q;*Wp — > to + iT) for momentum transfer Q = k/ — and energy transfer u>, where kj and k^ are the inci- 
dent and final neutron wave vectors, respectively. In high-Tc superconductors, such as the cuprates compounds, the 
phenomenological parameter T is a positive function of energy, momentum transfer, temperature, and doping, and 
could be measured by the width of the resonance. Since the position of the peak does not depend on T but does 
depend on the doping, we shall use the available ARPES data for underdoped Y Ba2Cu^,OQ.'j , assuming that T — > + . 
In this case the imaginary part of Xij is a delta function centered at the pole of the real part of the generalized spin 
susceptibility. Note that when T is a positive function, the denominator in the spin susceptibility vanishes if its real 
and imaginary parts vanish simultaneously. 

The generalized spin susceptibility, related to the Hamiltonian (J]), is defined as follows^: 



where 



Xij(Q;^ P ) = 2 / e l ^ u < T u {Mi(Q,u)Mi(-Q,0)} > du, 
Jo 



M,(Q, U ) = I {a t ) aal £ q e fi "^ +q/2 ,>Q_ qA .< e -^ 



(2) 



Here ipq,a = (1/N) £\ e l<l ' Ti ipi :a , the variable u ranges from to j3 = l/fc^T (T and fee are the temperature and 
the Boltzmann constant), and uj p — (2ir/j3)p;p = 0, ±1, ±2, .... From Eq. @ follows that the generalized spin 
susceptibility and the two-particle Green's function share common poles. 

The commensurate and incommensurate peak structures associated with the neutron cross section in YBaCuO 
have been studied within the Fermi-liquid-based scenario using the single-band Hubbard t — U model^— or the t — J 
modeliS^— . The techniques that have been used are based on (i) the Monte Carlo numerical calculations-, (ii) the 
random phase approximation (RPA) for the magnetic susceptibility^, (iii) the mean-field approximatio n 10 ' 11 and 
(iv) the RPA combined with the slave-boson mean field scheme^—. In the RPA the two-particle Green's function is 
replaced by a product of two single-particle zero-temperature Green's functions^ to obtain 

X (Q; W ) = X (0) (Q;w)/(l + U qX ^(Q;uj)), (3) 

where the BCS susceptibility x^°^(Q;cj) is given by the usual expression (in our notations x^HQ;^) = ^77) 



frv,A- 1 \^ri £ k£k+Q + A k A k+Q 
(Q >W )--2, J- E{k)E{k + Q) 



k 



E(k) +E{k + Q) 



(uj + iT) 2 - [£(k) + £(k + Q)] 2 



(4) 



and Uq — U for the Hubbard model, and Uq — —2J(cosQ x + cosQ y ) for the t — J model. Here £k is the mean-field 



electron energy , A k is the gap function and E(k) = y£k + A k . 

Since there is a consensus that the calculations based upon equation overestimate spin fluctuations because the 
RPA neglects the mixing between the spin channel and other channels. The coupling of the spin and two tt channels 
(a three-channel response-function theory) leads in the generalized random phase approximation (GRPA) to a set of 
three coupled equations. When the extended spin channel is added to the previous three channels, we have a set of 
four coupled equations (a four-channel theory)^. In what follows, the energy of the resonances are obtained from the 
solution of 20 coupled Bethe-Salpeter (BS) equations for the collective modes in GRPA, i.e. the resonances emerge 
due to the mixing between the spin channel and other 19 channels. 



II. BETHE-SALPETER EQUATIONS FOR THE COLLECTIVE MODES 

The interaction part of the t — U — V — J Hamiltonian is quartic in the Grassmann fermion fields so the functional 
integrals cannot be evaluated exactly. However, we can transform the quartic terms to a quadratic form by applying 
the Hubbard-Stratonovich transformation for the electron operators^: 



/ £Aexp i^ Q (z)^°j- 1 (z,z')A (9 (z) + ^( 2/ )r^(y;x|z)^)A Q (z) 



(o), 



exp 



-i^(y)fi 0) (y; x\z)^x)D^> (z, z')^{y')f^ (y'; 



)(0) 
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The t — U — V — J model has an additional term, the spin-dependent interaction Heisenberg interaction 
JJ2<i,j> Sj.Sj = Ji + J 2 , which consists of two terms: J\ = 7E<,, 3 >["i,t%t + ^i.4-%4- _ "i,t"j,4 _ ^»,4-%t] 

This requires to introduce four-component Nambu fermion 



and J 2 = |E<i,j> ^t^U^L^t + ^1,^,^1^3,1 



fields 



and 



ip(x) = (^p^(x)ilj^(x)^ (x)ipi(x)j , where x and y are composite variables and the field operators obey anticom- 
mutation relations. The Fourier transforms of the 4x4 matrices D^l and (a, ft = 1, 2, 3, 4) can be written in terms 

of the Pauli a u Dirac 7 and alpha matrices^: £> (0 ) = ( D ^ J V f ^ = (7 ±a z )/2 and f|°] = (a x ±ia y )/2, where 



<Ti 

OyCJiOy 



, L>i = (J(k) - F(k)) ctq + (£/ - - V(k)) a x and L> 2 = 2 J(k)a x . For a square lattice and 



nearest-neighbor interactions V(k) = 4V(cosk x + cosk y ) and J(k) = J(cosk x +cosk y ). Now, we can establish a one- 
to-one correspondence between the system under consideration and a system which consists of a four-component boson 

field A a (z) interacting with fermion fields rp(y) and ip{x). The action of the model system is S — S^f 1 + Sq A ^ + S^ - ^ 

where: S { e) = %)G^- 1 (y; x)$(x), = \A a {z)Dff\z,z')A p {z') and S^~ A ) ^%)W\y,x \ z)$(x)A a (z). 

Here, we have used composite variables x,y,z = {r^tt}, where is a lattice site vector, and variable u range from 
to f3 = l/ksT (T and ks are the temperature and the Boltzmann constant). We set H = 1 and we use the 
summation-integration convention: that repeated variables are summed up or integrated over. 

In Ref. [l9(, the spectrum of the collective excitations of the extended Hubbard model has been obtained from 
the Dyson equation for the boson Green's function D in terms of the proper self-energy. In what follows we shall 
obtain the spectrum of the collective excitations directly from the solutions of the BS equations for the two-particle 
Green's function. It can be shown that for a singlet superconductivity and d-wave pairing J\ and J2 terms contribute 
separately to the collective modes. Since the RPA expression ([3]) for the spin response function can be obtained by 
keeping only the J\ interaction term in our BS equations ([5]) and ([5]) when U=V=0, we shall neglect the contributions 
to the BS equations due to the J2 term. Following the same steps as in Refs. 20.21], we can derive a set of two BS 
equations for the collective mode w(Q) and corresponding BS amplitudes: 

KQ) - e(k, Q)]G+(k, Q) = £ E q [ 7 k, Q 7q,Q + k, Q Z q , Q ] G+(q, Q) 

-|Eq [7k,Q7q,Q - ? k,Q' q ,Q] G~ (q, Q) 

-27v E q [V(k - q) + J ( k - <l)] [7k,Q7 q ,Q + Iu,qI«,q] G+(q, Q) 
~27? E q I^( k - q) - ^( k - q)] [7k,Q7q,Q + m k , Q m q , Q ] G+(q, Q) 

+ 37? E q [V(k - q) + J(k - q)] [7k, Q7q ,Q - Zk.Q^q] G~(q, Q) 
+ 27? E q I^( k - q) - J ( k - q)] [7k,Q 7q ,Q - m k)Q m q)Q ] G"(q, Q) 
E q 7k,Q7q,Q (G + (q, Q) - G-(q, Q)) 
+ U ~l V N iQ) E q m k , Q m q>Q [G+(q, Q) + G~ (q, Q)] , (5) 

[w(Q) + s(k, Q)]G-(k, Q) = E q [7k,Q 7q ,Q + /k, Q ^ Q ] G-(q, Q) 
+ m E q [7k,Q7q,Q ~ ik,Q/q,q] G+(q, Q) 
+ 27? E q ^( k ~ q) + ^( k - q)] [7k, Q7q , Q + Ikq^q] G-(q, Q) 
+ 27? S q I^( k - q) - J ( k - q)] [7k,Q 7q ,Q + m k)Q m q)Q ] G"(q, Q) 

-27? E q m k - q) + ^( k - q)] [7k,Q7q,Q - ik,Q/q,Q] G+(q, Q) 
~27? E q I^( k - q) - J ( k - q)] [7k,Q7 q ,Q - m k , Q m q , Q ] G+(q, Q) 



U-2J(Q) 
2N 



E q 7k,Q7q,Q (G+(q, Q) - G-(q, Q)) 



- U ~l V N iQ) E q mk, Q m q ,Q [G+ (q, Q) + G~ (q, Q)] . (6) 

Heree(k,Q) — E(k + Q) + E(k), and we use the same form factors as in Ref. [2C| : 7 kj q = u k u k+ q +w k w k +Q, £ k .Q = 
WkMk+Q-fkWk+Q, 7k,Q = «k"k+Q-«k+q"k, and m kjQ = u k i; k +Q+w k +Qt> k where = = [1 +e(k)/S(k)] /2. 
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The Fourier transforms of V and J interactions are separable, i.e. V(k — q) = 2VXu^£ and J(k — q) 



JA k AT/2, 



and therefore, Eqs. ([5]) and ([6]) can be solved analytically. Here Ak = (sk, <ik, ssk, s^k) is an 1 x 4 matrix, and 
we have used the following notations: Sk = cos(k x ) + cos(fc y ), <7k = cos(k x ) — cos(k y ), ssk = sin(fc x ) + sin(fcj,) and 
ccik = sin(k x ) — sin(fc y ). Thus, the BS equations for the collective modes can be reduced to a set of 20 coupled linear 
homogeneous equations. The existence of a non-trivial solution requires that the secular determinant rfe<||x 1 — V\\ 

P Q 
Q T R 

V = diag{U, U, -(U - 2J(Q)), U - 2V(Q), -{2V + J/2), -(2V + J/2), -{2V - J/2), -{2V - J/2)) are 20 x 20 
matrices. Here, P and Q are 4x4 and 4 x 16 blocks, respectively, while R is 16 x 16 block (in what follows 
i,j = 1,2,3,4): 



is equal to zero, where the bare mean-field-quasiparticle response function x — 



and the interaction 



P 



I, 



7 7j7 



7,7 ^7^ 7,7 

Jl,J 



7 7j7 



J7,m 
J-7,m 



J~{,m Il,m Jy,m -^rr 



7? 



7:7 

7,' 



77 

7,m Z,m 



7* J ~ 

7,7 



in 

u 7,m 
1 Lm 



til- 
7:7 

VJ T l 3 

7,m m,m 



7~ J 

7,m 



p ~ 

7,7 



J 



J 7,7 J 7,' 

7,1 Ijl 

P - 7?_ ft _ 

7,7 i,7 7,7 

p ri 

7,m Z,m 



!>7 



J 



7~ 

7,m 



7,m 



7,m 



The quantities 7 ajfc = F Q:6 (e(k, Q)) and J Q:6 = F Q , 6 (w), the 1x4 matrices P a b = F£ b (e(k,Q)) and J* b = 7^ b (w), 
and the 4x4 matrices 7* J 6 = F^ b (e(k, Q)) and J* J 6 = F^ b [w) are defined as follows (the quantities a(k, Q) and 
Q) = ^k,Q, m k ,Q, 7k,Q or 7k,C}): 



7^) = -^ 



TV ^ uj 2 

k 



za(k,Q)6(k,Q) 
(k,Q) 



TV ^ 

k 



xa(k,Q)&(k,Q)A k 
2 -e 2 (k,Q) ' 



W = 77 



1 ^ za(k, Q)6(k, Q) ^ T 



A k A k 



III. NUMERICAL CALCULATIONS 



The sum over k can be replaced by a double integral over k x and k y , both from the first Brillouin zone. After 
that, we applied the substitutions x = tanfc^/4 and y — tanfc^/4 to rewrite the integrals in the form of Gaussian 
quadrature J —1 dx J _ 1 dyf(x, y)/(l + x)(l + y). The corresponding integrals are numerically evaluated using 49 x 49 

(xi,yj) points: dx /_ x dyf(x, y)/(l + x)(l + y) = Y,t=i Ei=i w i w jf( x U Vj), where u; 4 is the corresponding weight. 
The mean-field electron energy £k has a tight-binding form 



£k = — %t (cos fca; + cos k y ) + At' cos k x cos /% — 2i"(cos 2k x + cos 2fc y ) 



A 4 



(7) 



obtained by fitting the ARPES data with a chemical potential fi and hopping amplitudes ti for first to third nearest 
neighbors on a square lattice. Using the established approximate parabolic relationship 22 T c /T cjnax = 1 — 82.6(p — 
0.16) 2 , where Tc,max ~ 93K is the maximum transition temperature of the system, T c — 67K is the transition 
temperature for underdoped YBa2CusOe.7, we find that the hole doping is p = 0.10. At that level of doping the 
ARPES parameters are obtained in Ref^: t = 0.25 eV, t' = QAt, t" = 0.0444i and fj, = -0.27 eV. In the case 
of d-pairing the gap function is Ak = Atfk/2, where the gap maximum A should agree with ARPES experiments. 
In the case of underdoped Y Ba-iCu^O^j the gap maximum has to be between the corresponding A = 66 meV in 
Y Ba 2 CuzOQ,Q and A = 50 meV in Y ^a^CuaOg.g^, so we set A = 60 meV. The BCS gap equation is 



V,,, 



dk 



d' 2 

"k 



.(2tt) 2 E(kY 



where V^p = 2V + 3 J/2. The numerical solution of the gap equation provides V$ — 265 meV. 
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FIG. 1: Set of points in U, J parameter space which reproduces the commensurate resonance at 40 meV at point Q = (n/a, n/a) 
. The set is fitted with the linear formula U = -3.985 J + 1.01 (U and J are in eV). Note that V = 7^/2-3.7/4 where = 265 
meV is calculated from t he g ap equation by using the set of parameters given in Ref. [23|. The maximum value of the energy 
gap is A = 60 meV Ref. [Hf. 

During the second step of our approach we solved numerically the BS equations to obtain the spectrum of the 
collective modes w(Q) at the commensurate point Q = (n/a, it /a), as well as at four incommensurate points 7r/a(l ± 
S, 1) and 7r/o(l, 1 ± 5). Here 6 is the deviation of the collective mode position from Q . The collective-mode energy 
is the same at any of the four incommensurate points. Since the spin response function and the two-particle Green's 
function share the same poles, we can use the 40 meV solution of the BS equations at Q to obtain a relation between U 
and J parameters, which is presented in FIG. 1. It can be seen that the linear formula U = —3.985 J+ 1.01 agrees very 
well with the numerical results (U and J are in eV). By means of the last relation we have solved the BS equations for 
the one of the four incommensurate 24 meV peaks observed at (Q x ,Q y ) =7r/a(l±5,l) and (Q x ,Q y ) — vr/a(l, 1 ± 5), 
where S = 0.22i. The solution provides the following interaction strengths: J ~ 129 meV, V ~ 35.7 meV and U ~ 495 
meV. 

We have tested the above values of the interaction strengths by calculating the positions of the incommensurate 
peaks at 32 meV. The Bethe-Salpeter equations with the above strengths provide the deviation from Q of about 
5 = 0.196. This result is in agreement with the experimentally obtained deviation of 6 ~ 0.19 (see FIG. 2 in Ref. pj). 
The fact that our approach is able to reproduce the 40 meV peak as well as the two peaks at 24 meV and 32 meV 
allows us to conclude that the commensurate resonance and incommensurate peaks have a common origin. 

IV. SUMMARY AND DISCUTION 

The strength for J should be comparable to the strength of the superexchange interactions in the underdoped 
antiferromagnetic insulator state of the cuprates. The superexchange interaction in cuprates has been studied by 
using several experimental tools, and is now known to be not strongly dependent on materials with the magnitude of 
0.1 — 0.12 eV. Our value of J = 0.129 eV is larger than the accepted experimental values though some theoretical papers 
have predicted different magnitudes. For example, in Ref. [25| the value of J = 0.22 eV has been predicted using 
standard cuprate parameters. The calculated value of the superexchange interaction in Ref. }26i] is about J = 0.16 
eV. The differences could be due to the fact that calculations are very sensitive to the values of the superconducting 
gap and the tight-binding parameters, and therefore, somewhat different tight-binding form of the mean-field electron 
energy could bring the calculated J closer to the experimental value. 

In summary, we have demonstrated that the strengths of the interactions in cuprates can be obtained if we have 
the angle-resolved photoemission and inelastic neutron scattering data collected on the same crystals of the high- 
temperature superconductor. We do not wish to repeat the theoretical arguments that were advanced against the 
stripes model, but our unified description of the peaks based on the t — U — V — J model strongly supports the 
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hypothesis that the commensurate resonance and the incommensurate peaks in cuprates compounds have a common 
origin. 
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